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The temperature dependence of the nuclear matter isovector symmetry energy coefficient (Ao,i) 
is investigated in the framework of the generalized nuclear polarizability with Skyrme interactions, 
as worked out in Refs. P, @|- The variation of Aq,i(T) is very small (of the order of 1 MeV) 
for temperatures (T) in the range of and 18 MeV. Different behaviors with temperature are 
found strongly depending on the Skyrme parameterization, in particular at densities lower than the 
saturation density po. 
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I. INTRODUCTION 



The neutron-proton nuclear symmetry energy is cur- 
rently under intense theoretical and experimental inves- 
tigations not only due to its relevance to different as- 
pects of the nuclear structure and dynamics but also to 
description of dense stars structure and of the Super- 
novae mechanism 0, i, & H H 0, S, IU • Actually, the 
possibility of extracting better experimental values for it 
at different densities and temperatures is under constant 
improvement. In this respect, (multi)fragmentation pro- 
cesses in heavy ion collisions provide experimental knowl- 
edge about the liquid-gas phase transition and therefore 
about the nuclear forces at low nuclear densities and 
excitation energies. For example, the isoscaling found 
in yields of multi-fragmentation experiments 0] depends 
directly on the isovector symmetry energy (coefficient) 
which was found to decreases considerably with the ex- 
citation energy according to experimental data [ll[. In 
fact, in experimental conditions it is very difficult to ex- 
tract unambiguous behavior with each of the thermody- 
namic observables involved (T, p) [Tl| and in earlier de- 
scriptions of experimental results some groups have con- 
sidered a seemingly too strong variation of the symmetry 
energy with temperature In different (more recent) 
analysis it has been found that the strong decrease of 
the symmetry energy in experimental conditions should 
also be due to to the excitation energy dependence of 
surface effects [l5[ and to the expansion of the system 
which implies the lowering the total density 0, [lyj . Be- 
sides that, recently De and Samaddar [ijj have argued 
that the symmetry free energy (l3j is the parameter that 
appears in the scaling of multi-fragmentation. However 
their arguments stand for the specific analysis of multi- 
fragmentation processes and we intend to discuss rather 
the behavior the symmetry energy coefficient behavior. 
Concerning the specific dependence on the temperature 
(up to T ~ 20MeV) there are several theoretical indica- 



tions of very weak decrease [13, [17|, [18|, [19|, l2fj although in 
some works, for finite nuclei and nuclear matter, a small 
(or very small) increase was found [l], HM HH • In some of 
these works no unique tendency was found mainly for dif- 
ferent finite nuclei, for example in Ref. 18]. In Ref.[l] few 
preliminary numerical results showed a very small vari- 
ation of the nuclear matter symmetry energy coefficient 
with temperature using Skyrme forces. In the present 
work we perform a quite extensive investigation of this 
subject with the generalized polarizability as proposed in 
Refs. [3, ||. 

In usual mass formulae the binding energy per nucleon 
depends on the n-p asymmetry with the following form: 
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where N, Z, A are the neutron, proton and mass num- 
bers, the isovector symmetry energy coefficient (s.e.c.) is 
a T . This coefficient is a measure of the energy needed 
to increase n-p asymmetry. Different powers of the n-p 
asymmetry (N — Zy /A (J ^ 2) are quite smaller, al- 
though the n-p asymmetry modifies many other terms 
22]. The value of the n-p s.e.c. in fits of the mass formu- 
lae for stable nuclei is quite well known Q. For an infi- 
nite medium, roughly speaking, nucleon numbers (N, Z) 
might be replaced by the corresponding densities p n ,p p 
for many purposes depending on the volume occupied by 
each nucleon species, see for example Ref. |23{ . 

A particular interesting way of obtaining the symme- 
try energy coefficient, eventually in different conditions, 
was found by the nuclear isovector polarizability Aq i 
[U [Hj]. Generalized nuclear polarizabilities and their 
dependencies on several parameters simultaneously (such 
as total density p, exchanged momentum and energy q, ui, 
neutron-proton asymmetry b = ^ — 1) were investigated 

quite extensively with Skyrme-type forces @, [H, [2(| H3] • 
Deviations from the quadratic form of expression ((T|) 
might appear depending on the particular density fluc- 
tuations considered for a given set of nucleon densities 
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for these calculations with Skyrme density functionals. 
Recent investigations have revealed further relations of 
these (non relativistic) effective interactions with rela- 
tivistic models (HIH- 

This work exhibits the temperature dependence of the 
isovector generalized s.e.c, as proposed in Refs. [HQ us- 
ing Skyrme type forces. For that, the generalized isovec- 
tor screening function with different parameterizations of 
Skyrme effective interactions are shown at different den- 
sities and n-p asymmetries. 



II. 



GENERALIZED POLARIZABILITIES 



For an asymmetric medium in n-p densities, the gen- 
eralized screening function using Skyrme forces, for zero 
energy and momentum exchange (uj = q = 0), can be 
written in a compact notation, as [HQ: 



A s ,t = ^ {l + 2vt t] N c + wt' t] M;{ Pc + Pd )+ 



+ 12M;V 1 {s ' t) V (s ' t) (N cPd - p c N d ) + (2) 
+ (vf ' t] f (m(M;fp cPd - l6M*M c N d ) | . 

Where Vo* and are functions of the Skyrme 

forces parameters in a given channel of the effective 
nuclear interaction of (spin, isospin) denoted by super- 
scripts (s,t). Therefore they carry the main contribu- 
tions of the effective NN interaction, distinguishing each 
of the channel of the particle-hole interaction. In partic- 
ular, the function V^ ,l) = (i 2 (l + 2as 2 ) -ti(l + 2jci))/16 
is a combination of (momentum dependent) Skyrme pa- 
rameters that contributes to the usual nucleon effective 
masses in the framework of Skyrme calculation, m* np . 
The effective masses of neutron and protons are func- 
tions of the total density and of the neutron and proton 
densities, see for example in (24[. On the other hand, 
the function Vjj depends mainly on to and £3. While 
yt ' 1 ) of each of the Skyrme parameterization used in this 
work do not have meaningful different values, the func- 
tion Vi has considerably different values because it 

might be zero, i.e. for some forces V - / ' 13 = 0. Although 
relevant these differences will be shown to not be enough 
to provide too much different behaviors for A s ,t(T). The 
mixed functions p v , M v and N v reduce to densities of the 
Skyrme-Hartree-Fock approach at zero temperature, be- 
ing respectively the nucleon and kinetic energy densities 
and the densities of states. They are the zero frequency 
and zero momentum generalized Lindhard functions cal- 
culated in Ref. [l[. The total densities, p,N, are writ- 
ten without any index. These functions (p v , M v , N v ) are 
given respectively by: 



Pv = V Pn + (1 - V)p p , 

M v = vM n + (1 - v)M p , 
N v = vN n + (1 - v)N p . 



(3) 



In these expressions v stands for two different n-p asym- 
metry coefficients (c, d) , being that we made use of the 
following asymmetry coefficients: 
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Therefore, by fixing the parameter of density asymmetry, 
6 (for a given total nuclear density), the neutron and 
proton densities are found as well as the other asymmetry 
parameters (a, c, d) for a given Skyrme force [H, 

There is a further mass parameter in expression 
M* which is in fact a kind of reduced mass. It is given 
by 
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Although the calculation has been carried out in such 
a way to provide final expressions for each of the chan- 
nels of the particle-hole interaction (isovector, spin, spin- 
isovector and scalar) only the n-p one will be investigated 
in this paper. 



A. Varying Temperature 



The densities N a ,p a and M a (for a = n,p neutrons 
and protons) are the basic input for the temperature de- 
pendence of the polarizabilities. At finite temperature 
these functions are given respectively by integrals writ- 
ten as: 



(JV ai 3p Q , 4M a ) 
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df a (k)(k.m* a , k A - fc 5 ). (4) 



In these expressions, df a (k) is the measure of integration 
in terms of the usual free fermion occupation numbers 
fa(k) for neutrons and protons (a = n,p). At T = for 
the usual Fermi occupation number we have: 



df a (k) = -6(k 



(5) 



Where k^ is the Fermi momentum for each of the nu- 
cleon species. In this case the integration is trivial. 
Therefore all the temperature dependence of the density- 
like quantities p, N, M is encapsulated in the integrals 
above Q which can, at most, yield smoother results for 
the integral. Furthermore N V (T), M V (T) are the only 
parameters that vary with temperature since p a are kept 
constant. The zero temperature limit was considered pre- 
viously 0, i] and it shows more explicitly, as mentioned 
above, the effect of each of the Skyrme force parameters 

(s t) { s t) 

through the functions Vq ' and V+ ' as well as the ef- 
fective masses. This issue is extremely relevant for the 
resulting Ao,i(T). 

From the general expression @ an useful (simplified) 
limit is recovered in which the behavior with tempera- 
ture can be understood in detail. For instance, consider 
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the limit in which the function V±° appears only in 
the leading order of the symmetric n-p function. This 
is achieved with b = a = and d = c = 1/2, yielding 
to* = m* = to* = M* and p p = p n = p/2. We obtain 
an expression of the following form, for (s,i) = (0, 1), in 
the n-p symmetric limit: 



± + 2Vr ) +Grn*Vr%)+h.o.(Z) 



Where h.o. stands for the higher order terms in T^- ' 1 ^. 
This expression reproduces exactly a T , which is the usual 
symmetry ener gy c oefficient in the Skyrme-Hartree-Fock 
approximation [24|. The first and third terms of a T de- 
pend on N and therefore they show a very small variation 
with temperature only due to N(T) according to expres- 
sion (H]). The behavior of the function N a (T) is mono- 
tonic and it decreases with T. This is the main feature 
for understanding the numerical results of expression © . 

The higher order terms, in the limit of n-p symmetric 
matter, are given by: 



h.t 



(1/ 1 (0 ' 1 ) 2 (9(to*) 2 P 2 -8to*M7V)-^. (7) 
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In these terms, and mainly for n-p asymmetric matter (in 
the complete expression @), the imbalance between the 
T-dependence of M v (T) and N v (T) determines whether 
the increasing behavior of N(T) with temperature is the 
leading one or not. Basically this is seen from the (over- 
all) denominator of expression ((21) by reminding that 
N(T) is a decreasing function of the temperature in the 
Skyrme-Hartree-Fock level, and on the other hand M V (T) 
is an increasing function of the temperature. The differ- 
ent resulting behaviors might appear due to the relative 
values of their coefficients, i.e. Vq ^ and V^ ' . This 
is noticed in the results exhibited in the next section. 
Although it might be expected that for asymmetric n-p 
matter the variation of the polarizability is larger since 
there are more terms at work, this will be shown to be 
not really sizeable for the forces considered in this work. 
The complete expression is quite complicated such that it 
might not exhibit a simple and unique behavior in more 
general situations. 

It is worth to mention some recent results claiming 
that the symmetry free energy is the quantity that re- 
ally rules multi-fragmentation pj|, [H, H(J. Basically 
this corresponds to considering the entropic contribu- 
tion which amounts basically to extra additive terms. 
This can be qualitatively seen as follows. The calcu- 
lation of the polarizability in such case should depart 
from a free energy in the presence of a n-p asymmetry 
and of an infinitesimal external source (e) that induces 
fluctuations of the densities of neutrons from protons, 
i.e: T{p + 5p np ) = £(p + Sp np ) -TS(p + Sp np ) + eSp np , 
where Sp np = Sp n — 5p p . The entropy can be expanded 
in terms of p n — p p to make explicite its contribution to 
the symmetry (free) energy: F(T) — £ (T) + a T (5p np ) 2 + 
SW(T)6p np - S^{T)(Sp,[ p ) 2 +eS Pnp + .. Where S«(r) 



are the leading contributions of the entropy for the sym- 
metry free energy. As it was shown in Ref . Q , the linear 
term in 8p np might be incorporated into the usual cal- 
culation (considering only the quadratic terms) and it 
is not considered explicitly below. The polarizability is 
then given by: 
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The final symmetry free energy coefficient can be written 
as a{{T) = a T (T) - S< 2 )(T).' Therefore we can expect 
that the entropic contribution would appear mainly as 
additive terms for the screenin g fu nction. This is seen 
in the results of Refs. pj|, [l3l l3C)j|. However a micro- 
scopic investigation of this quantity, with its eventual 
relevance for the multi-fragmentation processes, is out- 
side the scope of the present work. 



B. Results 

The temperature dependence of the isovector polariz- 
ability Aq.i(T) is shown in figures 1 to 4 for the following 
Skyrme forces: SGII from Ref. [H[, SLyb from Ref. [if 
(which is sometimes referred to as SLy4 in the litera- 
ture )and two parameterizations SkCS4 and SkSC6 from 
Ref. [33| • These last two Skyrme parameterizations have a 
slightly more intricate density dependence although the 
resulting functions Vi° are zero. The zero function 



(0,1) 



(that carries the main part of the momentum 



dependent Skyrme forces) brings a lot of simplification 
in the dependence on the temperature as discussed in 
the last section and it is noticed in the figures below. In 
these cases the behavior of A s ,t(T) is always monotonic 



(0,1) 



Differ- 



depending on N(T), and slightly less on V c 
ent total densities and n-p density asymmetries are also 
considered. 

In Figure 1 the function Aq,i(T) is shown at the sat- 
uration density, p = po, and zero n-p density asymmetry 
(b = 0) with the following Skyrme parameterizations: 
SkSC4 (circles) SkSC6 (x) SGII (squares) 

[3ll ] . SLyb (diamonds) [32l |. The variation with tem- 
perature (up to T ~ 18 MeV) is quite small, reaching 
A.Ao,i — 0.5 — 1.0 MeV depending on the interaction, 
and even nearly zero for the SGII force. The slope is al- 
ways positive, although smaller at high temperatures. As 
noticed after expression ^ the function N(T) decreases 
with temperature within the Skyrme-Hartree-Fock ap- 
proach. The force SLyb is the one with larger variation 
in Anj(T). This trend of small variation was found be- 
fore pi, although it disagrees with the small decrease of 
the s.e.c. with T found in different works for lower den- 
sities [H, [U H [H IM EH- However it is also worth 
to point out that in these references the s.e.c. was in- 
vestigated in the regime of very low total density. The 
variation found for the range of T — up to T ~ 15MeV 
is not large in all these works, and it can be of the order 
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of 1 up to 3 MeV (though negative), a little bit larger 
than the present results. 

Some further remarks to understand the behavior with 
temperature are in order. The chemical potential fixes 
the nucleonic density which is kept constant for all tem- 
peratures. In the calculations with Skyrme interactions, 
the kinetic part of a T and the terms with N(T), M(T) in 
expression (j2J|) are temperature-dependent. In particular 
whereas the functions N a slightly decrease with the tem- 
perature, the densities M a increase slightly. Having this 
in mind and analyzing expression {2} we can expect that 
n-p asymmetry (6, d, c ^ 0) favors different behaviors of 
«4o,i(T). The relative variation of the potential energy 
part of the symmetry energy at the Hartree Fock level is 
nearly zero since it depends mostly on p at not very low 
densities. 

The same Skyrme parameterizations (and symbols) are 
used in Figure 2. In this figure, Aq^{T) is exhibited for 
p = 0.75po with b — 0.25 (full symbols and x) and 6 = 
(empty symbols and +). For these cases (of lower nuclear 
matter density), the behavior of the polarizability with 
temperature is non-monotonic in the case of force SGII 
(squares). For Skyrme forces SkSC4, SkSC6 and SLy, 
the polarizability Aq,i(T) is nearly constant at very low 
temperatures and it starts increasing smoothly around 
T ~ 5MeV until T ~ 10 - 15MeV. For the force SGII, 
Ao,i decreases for low temperatures and smoothly in- 
creases for temperatures higher than nearly 3 — 5McV. 
By comparing the relative variation of the polarizability 
for n-p symmetric and for b = 0.25 n-p asymmetric mat- 
ter, we find no further meaningful difference. It is worth 
to emphasize that the lower total density makes possible 
this non monotonic behavior. This is produced the differ- 
ent behaviors with temperature of the functions N a (T) 
and M a (T) for given V^ 0,1] and v/ 0,1) . This becomes 
clearer with Figure 3, where p — 0.5po- As noticed above, 
the n-p asymmetry might amplify this non-monotonic be- 
havior although the difference is very small in the cases 
we show. 

The same Skyrme parameterizations (and symbols) are 
used in Fig. 3, where Ao t i(T) is exhibited for the still 
lower density p = 0.5po with: b = 0.5 (full solid symbols 
and x) and b = 0.25 (empty symbols and +). Differently 
from all the results shown above the only Skyrme force 
that exhibit the non-monotonic behavior is SLyb for both 
b = 0.5 and b = 0.25. Furthermore we notice that the 
stronger variation for all the forces (even if they are very 
small of the order of .5MeV) occur below T ~ 5 MeV 
or T ~ lOMeV. Above these temperatures, the isovector 
polarizability variation is smaller. The reason why the 
SGII and SLyb forces have the non-monotonic behavior 
at different densities is explained by the different values 
of the functions Vq ' 1 ^, in terms of the to, t% Skyrme pa- 
rameters 0, 0, and also v} ' 1 ^ (non zero). 

For the range of lower nuclear densities some fur- 
ther conclusions can be extracted by comparing the fig- 
ures 2 and 3, in particular for b = .25, (which means 



p n = 1.25p p ). We notice that the behavior of Aqa{T) 
is different depending on the effective force. While the 
parameterizations SkSC 4 and SkSC 6 do not provide 
any different behavior (apart from an eventual overall 
total variation of Aq.i(T)), the forces SGII and SLyb 
present different trends for p = .5po (Fig. 3: empty 
squares, SGII, and empty diamonds, SLyb, respectively) 
and p = .75 po (Fig. 2: full squares and full diamonds re- 
spectively). The non-monotonic behavior appears for the 
SGII parameterization at p — .75po while for the SLyb 
one it appears when p = .5po. As discussed in the be- 
ginning of the last section, because of the complicated 
form of the expression ([2]) and of the behavior of the 
functions N(T) and M(T), the results from the polariz- 
abilitics with Skyrme forces are not always monotonic. A 
suitable quantity for comparing the results from different 
forces, and even different methods, is the total variation 
AAo.i = A ,i(T = 20) - A ,x(T = 0). This quantity 
is (quite) small in all the works of the field, with small 
differences also due to the particular nuclear density un- 
der consideration being also seen also in other works us- 
ing different approaches 0, [l7l [3 H3| as analyzed in 
|26| . One of the main outcomes of these comparisons goes 
along with the above remarks: the behavior of Ao,i{T) 
might be different at po, 0.75po and 0.5po, eventually for 
still lower densities analyzed in other works. 

In Figure 4, Ao^(T) is exhibited for a density higher 
than po, i.e. p = 1.33po with b = (full symbols and x) 
and b = 0.25 (empty symbols and +). The same kind 
of behavior found for p = p$, in Figure 1, is present in 
Figure 4. The isovector polarizability very smoothly in- 
creases with temperature, although the variation is con- 
siderably smaller at high temperatures. 

This analysis suggests that the T-dcpcndcnce of the 
s.e.c. is strongly dependent on the nuclear matter den- 
sity p, and it is also suitable for shedding light on the 
nuclear effective interactions expected to be reliable with 
good predictive power. Nevertheless we emphasize that 
experimental data with temperature are very difficult to 
be extract unambiguously [111 ]. 



III. SUMMARY 

To summarize we conclude that the bulk isovector sym- 
metry energy does not vary considerably in a quite wide 
range of temperatures within the isovector polarizability 
with Skyrme forces. As noticed in other works the de- 
pendence of Aq.i on the n-p density asymmetry is prob- 
ably too strong. This issue, on the other hand, does not 
modify the variation with the temperature meaningfully. 
In the framework of the Skyrme-Hartree-Fock parame- 
terization, temperature effects arise from the functions 
M, N given by expressions Q . They depend on the gen- 
eral properties of Skyrme-Hartree-Fock approach. The 
variation of Ao t i(T) depends strongly on the particular 
Skyrme force, being always very smooth and small. The 
different contributions of the potential and kinetic parts 
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of the symmetry energy for fixed densities, by means of 
the functions N(T) and M(T), as well as the relative 

values of Vq and v} ' 1 ^, are responsible for these dif- 
ferent results of each Skyrme interaction. The variation 
(decrease) of N a (T) (which reduces to the n,p densities 
of states at zero temperature) with temperature is how- 
ever the most relevant contribution for the results. The 
larger variation of A .i(T) occurs for p < p Q , depend- 
ing strongly on the effective force parameterization (for 

which Vi (0,1) ^ 0, i.e. SGII and SLyb). The final be- 
havior is not always monotonic with temperature. The 
differences in the overall variation of Ao,i(T) for different 
densities below pojjs seen also in other works using dif- 
ferent approaches [Il,[l7],[ljl[Iil as pointed out and com- 
pared in Ref. [26j]. The trends exhibited by the Skyrme 
parameterizations suggest that the eventual experimen- 



tal knowledge of the behavior of the symmetry energy 
with the temperature will also contribute to better fine- 
tuning of the effective interaction as well as to improving 
its predictive power. This is clearer in the comparison 
between figures 2 and 3 (mainly for Skyrme forces SGII 
and SLyb). 
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FIG. 2: The inverse of isovector polarizability, »4o,i, as a 
function of the temperature for the same Skyrme forces as 
in figure 1. Considering: full symbols and symbol (x) for 
(p/Po = 0.75 and b — 0.25) and empty symbols and + for 
(p/po = 0.75 and 6 = 0.). 
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FIG. 1: The inverse of isovector polarizability, *4o,i, as a 
function of the temperature for different Skyrme forces: SkSC 
4 H3 (circles), SGII Q (squares), SLyb H2 (diamonds), 
SkSC 6 H (x or +). For: (p/po = 1 and b = 0). 
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FIG. 3: The inverse of isovector polarizability, *4o,i> as a 
function of the temperature for the same Skyrme forces as 
figure 1. Considering: full symbols and symbol x for (p/po = 
0.5 and b = 0.5) and empty symbols and + for (p/po = 0.5 
and b = 0.25). 
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FIG. 4: The inverse of isovector polarizability, Ao,i, as a 
function of the temperature for the same Skyrme forces as 
figure 1. Considering: full symbols and symbol x for (p/po = 
1.33 and b = 0.) and empty symbols and + for (p/po = 1-33 
and b = 0.25). 



